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Abstract. The average flux decrement shortward the Ly^ 
emission, due to the well-known "forest" of absorptions, has 
been measured in the spectra of 8 quasars. Quasi-simultaneous 
optical and lUE observations of the two low redshift quasars 
PKS 0637-75 (z=0.654) and MC 1104-^16 (z=0.632) have been 
carried out, obtaining relatively high S/N, spectrophotometri- 
cally calibrated data on their energy distribution from the rest 
frame H^ to the Lyman continuum. Six more quasars in the 
redshift range 2.5-3.4 have been observed in the optical do- 
main. For all the quasars the "intrinsic" continuum slope and 
normalization have been estimated longward the Ly^ emission 
and extrapolated towards the Lyman continuum to measure 
the average depressions, which have been compared with the 
model statistics of the Lyo, clouds. When all the known classes 
of absorbers are taken into account with plausible values for 
their equivalent width distribution and evolution, a good agree- 
ment is obtained with the observations. The results for the 
observed continuum decrement at 2 '^ 0.65 are identical to 
those predicted by the evolution with redshift of the number 
of Lya forest systems including the HST data and within 2a 
of the predicted value using the "standard" Ly^ evolution (as 
determined only at high z). 
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1. Introduction 

Considerable information about non-luminous extragalactic 
matter has been gained by studying the absorption features 
in quasar spectra. Particularly interesting in this respect is the 
issue of the evolution of the absorbers responsible for the "Ly^ 
forest systems". It is generally accepted that they are origi- 
nated by intergalactic primordial clouds (Sargent et al., 1980) 
and as such distinguished - although recent observations tend 
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to blur this difference - from the less numerous "metal-line 
systems" . 

Optical spectroscopy at high resolution available for a few 
QSOs at 2 ~ 3 shows a high number of hydrogen absorption 
lines whose statistical properties can be represented, in terms 
of the number density of hydrogen clouds, by a power-law in 
redshift with index ~ 2, indicating a cosmological evolution of 
the absorbers. 

At low resolution this "forest" of absorptions appears as 
a steepening of the continuum slope in the region shortward 
the Ly^ emission as compared to the slope defined on the red 
side of Lye (e.g. Steidel & Sargent, 1987). The average flux 
decrements Da, between Ly^ and Ly^ emissions, and Db, be- 
tween Ly/3 and the Lyman limit have been introduced to study 
the average properties of the absorption (Oke and Korycansky 
1982, Steidel and Sargent 1987). These quantities are, within 
relatively large limits, independent of the instrumental reso- 
lution, at least as far as the true continuum can be reliably 
extrapolated from the unabsorbed (red) region. 

Using the material available in the literature (Bechtold et 
al. , 1984; O' Brien et al. , 1988; Oke and Korycansky, 1982; 
Schneider, Schmidt and Gunn, 1989; Steidel and Sargent, 1987) 
together with new optical data obtained by us at high red- 
shifts, we have tried to disentangle the evolutionary details of 
the Ly^ absorbers (Giallongo & Cristiani, 1990). Excess ab- 
sorption with respect to the 7 ~ 2 evolution is found both at 
high and low redshifts (as could be expected since the evolu- 
tionary law has been normalized at z~ 3 where the bulk of 
the optical data has been obtained). The presence of differ- 
ent populations of absorbers might be suggested as well as a 
change in the physical conditions of the clouds. In any case 
the slower evolution found at low redshifts resembles that of 
metallic lines, perhaps indicating a link with the evolution of 
normal galaxies. 

The acquisition of new material is important especially at 
low redshifts, because in this region the uncertainties are still 
rather large and the constraints to the functional form of the 
evolution of the absorbers imposed by the data at the lower 
extreme of the observable redshift range can be particularly 
strong. For this reason we have observed 2 bright QSOs (V 



Fig. 1. lUE spectra of PKS 0637-75 and MC 1104+16, fluxes arc in units of 10" 
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~ 16) with z~ 0.65 (an especially convenient redshift for lUE 
observations since the QSO Ly^ emission falls then in the re- 
gion of overlap between the SWP and the LWP cameras) to 
provide additional good quality spectra needed to clarify the 
situation. In addition, optical observations of 6 quasars have 
been carried out in the optical domain, in order to improve the 
statistics and the quality of the data at z ~ 2.5 — 3. 



2. UV Observations 

Two consecutive lUE low-resolution long-wavelength spectra 
and one short-wavelength spectrum of PKS 0637-75 have been 
obtained on Dec. 25-26 1990 at VILSPA. One SWP and one 
LWP spectrum secured during a subsequent run on April 6- 
7 1991 allowed us a full coverage of the lUE spectral range 
(AA 1200+3300 A) also for MC 1104+167. Since the visual 
magnitude of both objects (V~16) is beyond the imaging ca- 



Table 1. Log of used lUE observations. 



Object 



Spectrum Obs. Date Exp. T. (m.) ECC 



0637-75 


SWP 10832 


80 Dec 12 


416 


— 


0637-75 


LWP 13013 


88 Apr 10 


120 


335 


0637-75 


LWP 13016 


88 Apr 11 


105 


335 


0637-75 


LWP 13022 


88 Apr 12 


120 


338 


0637-75 


LWP 19471 


90 Dec 25 


180 


331 


0637-75 


LWP 19472 


90 Dec 25 


215 


341 


0637-75 


SWP 40461 


90 Dec 26 


404 


331 


1104+167 


SWP 41309 


91 Apr 06 


400 


352 


1104+167 


LWP 20094 


91 Apr 07 


360 


452 



pabilities of the lUE FES, a standard blind offset technique 
from nearby stars was adopted to put the target into the slit. 



Fig. 2a. Optical spectra of Q0053-4029, Q0126-4050, Q0249+0222, Q0254-0137, relative fluxes are in units of 10" 



erg s -*- cm 



1-2 A-1 



The exposure times spanned from 3 to 7 hours giving a 
S/N ratio 10 ~ 15 in the best exposed regions. Taking into 
account that PKS 0637-75 did not show significant UV vari- 
ability over a decade, four pre-existing lower-quality spectra 
collected at GSFC between 1980 and 1988 have been included 
in our analysis. Their flux levels have been renormalized to 
the level observed in Dec 1990 and a weighted mean of all the 
data has then been carried out. A previous SWP spectrum 
of MC 1104+167 does not contain signal and was discarded. 
A complete log of used observations is reported in Table 1. 
This table also includes, when available, the VILSPA Expo- 
sure Classification Code (ECC) which gives an approximate 
measure of the image quality. 

Both our spectra and previous data have been consistently 
re-extracted from the so-called line~by~line spectrum provided 
by the standard lUESIPS processing. These spectra are already 
wavelength calibrated and photometrically linearized. Our sub- 
sequent extraction includes removal of cosmic ray hits and cam- 
era reseau marks as well as subtraction of a smoothed nearby 
background giving a time-integrated net spectrum. The final 
absolute flux calibrated spectrum is obtained by multiplying by 



the proper calibration function and dividing by the exposure 
time. Since our data have been processed with the current ver- 
sion of IUESIPS, the present LWP ITF2 calibration function 
(Cassatella et al. 1988) has been adopted for long-wavelength 
spectra. The short-wavelength data have been calibrated by 
means of the standard SWP function (Holm et al. 1982). A 
test extraction including the signal weighting on the basis of 
the observed PSF did not improve significantly our final S/N 
ratio. 

Both SWP and LWP data have been finally corrected for 
the time-dependent camera sensitivity degradation as pro- 
posed by Bohlin and Grillmair (1988) and Teays and Garhart 
(1990) respectively. The estimated absolute fiux error does not 
exceed 0.05 dex in both cameras. The lUE spectra of both 
quasars arc shown in Fig. 1. 



3. Optical Observations 

The journal of the optical spectroscopic observations carried 
out at the La Silla Observatory is given in Table 2. 



Table 2. Journal of the optical spectroscopic observations 



Name R.A. Dec. Date Telescope Spectrograph Slit width Resolution Sp. Range 

(1950.0) (1950.0) (arcscc) (A) (A) 



Q0053-4029 


00 53 51.5 


-40 29 30 


02 Sep 1987 


3.6m 


EFOSCl 


2.0 


15 


3750-7000 


Q0126-4050 


01 26 13.6 


-40 50 05 


03 Sep 1989 


2.2m 


B.&C. 


2.0 


20 


3400-8650 


Q0249+0222 


02 49 42.4 


+02 22 57 


30 Nov 1989 


3.6m 


EFOSCl 


2.0 


15 


3700-7000 


Q0254-0137 


02 54 07.9 


-01 37 49 


30 Nov 1989 


3.6m 


EFOSCl 


2.0 


15 


3700-7000 


Q0256-0000 


02 56 31.8 


-00 00 29 


02 Jan 1984 


3.6m 


B.&C. 


3.0 


20 


4350-8320 


Q0301-0035 


03 01 07.7 


-00 35 01 


17 Oct 1984 


3.6m 


B.&C. 


2.0 


20 


3700-8750 


PKS0637-75 


06 37 23.3 


-75 13 38 


20 Dec 1990 


3.6m 


EFOSCl 


1.5 


15 


3700-7000 


PKS0637-75 


06 37 23.3 


-75 13 38 


19 Dec 1986 


2.2m 


B.&C. 


6.0 


30 


3500-10000 


MC 1104+16 


11 04 36.7 


+16 44 17 


13 Apr 1991 


1.5m 


B.&C. 


2.0 


22 


4800-8800 



Fig. 2b. Optical spectra of Q0256-0000, Q0301-0035, PKS 0637-75, MC 1104+16 



For the reduction process the standard MIDAS facilities 
available at the Padova Department of Astronomy and at ESO 
Garching have been used. The raw data were sky-subtracted 



and corrected for pixel-to-pixel sensitivity variations by divi- 
sion by a suitably normalized exposure of the spectrum of an 
incandescent source. Wavelength calibration was carried out 



Fig. 3. log{v) - log{F^) combined spectra of PKS 0637-75 and MC1104+16. 



by comparison with exposures of Helium and Argon lamps. 
Absolute flux calibration was finally achieved by observations 
of standard stars listed by Oke (1974) and Stone (1977). 

For both PKS 0637-75 and MC 1104+16 the optical spec- 
trophotometry has been obtained quasi-simultaneously with 
respect to the lUE observations in order to minimize the ef- 
fects of variability. Three spectra of PKS 0637-75 obtained on 
December 1986 were normalized to the flux level of Dec 1990 
and merged with the 1990 spectrum. 

The optical spectra of all the quasars are shown in Fig. 
2a-b. 



4. Measuring the average flux decrement 

All the spectra have been corrected for galactic extinction ac- 
cording to Burstein & Heiles (1982) and rebinned to the QSO 
rest frame on the basis of their emission redshift, as shown in 
Fig. 3, in a log[v) — log{Fv) form. 



The slope of the continuum has been estimated in the form 
of a power-law fi, oc v~" (the dashed line in Fig. 3) on the 
basis of the "true continuum" regions at lower frequencies with 
respect to the Ly^ emission, and extrapolated to the regions A 
and B, as defined by Oke and Korycansky (1982). The average 
flux decrements Da and Db have been computed according to 
the equation: 



Di —< 1 — fobs/ fit 



>, 



i = A,B 



(1) 



where fobs and Unt are the observed and "intrinsic" fluxes per 
unit wavelength in the QSO rest frame. 

As emphasized in Giallongo and Cristiani (1990), the eval- 
uation of the flux decrements Da and Db is subject to many 
sources of uncertainty, closely related to the difficult estimate 
of the "intrinsic continuum" . At low redshifts these problems 
have relevant consequences on the interpretation of the data 
since, as shown below, the derived uncertainties (contrary to 



Table 3. Summary of the relevant data 



Name 


'^e.m 


Es-v 


dv 


Ba 


Da 


Db 


Bb 


i.s. Mg II EW 






(mag) 


(A > 1250 A) 


(1050-1170 A) 


(1040-1190 A) 


( 920-1015 A) 


( 912-1020 A) 


(A) 


PKS 0637-75 


0.654 ± 0.003 


0.12 


0.30 ± 0.05 


0.07 ±0.04 


0.06 ± 0.04 


0.18 ±0.05 


0.21 ±0.05 


6.5 ±0.5 


MC 1104+16 


0.632 ± 0.002 


0.02 


0.48 ± 0.05 


n 04 + "■''2 

U-U4 _ 0.01 


02 + ""2 
U-U^ _ 0.01 


11+ O''^ 
U.ii _ 0.02 


10 + O''^ 
U.iU _ 0.02 


4.0 ±0.3 


Q0053-4029 


2.758 ± 0.005 


0.01 


0.70 ± 0.08 


0.18 ±0.04 


0.18 ±0.04 








Q0126-4050 


2.594 ± 0.005 


0.01 


0.35 ±0.10 


0.23 ± 0.03 


0.19 ±0.03 








Q0249+0222 


2.805 ± 0.005 


0.06 


0.82 ±0.08 


0.17 ±0.03 


0.14 ±0.03 








Q0254-0137 


2.684 ± 0.005 


0.04 


0.53 ±0.05 


0.20 ± 0.02 


0.18 ±0.02 








Q0256-0000 


3.367 ± 0.004 


0.06 


0.61 ±0.05 


0.26 ± 0.02 


0.22 ±0.02 








Q030 1-0035 


3.205 ± 0.005 


0.06 


0.54 ±0.05 


0.30 ± 0.03 


0.27 ±0.03 











the situation at high redshift) are of the same magnitude as 
the measured effect. 

Composite quasar spectra (Cristiani and Vio 1990, Francis 
et al. 1991), used by us as a guide in the procedure of ex- 
trapolation, illustrate how few regions can be considered as 
indicators of the "true" continuum. In particular the excess 
emission around C III] (the blended Fe II 2000 A complex) is 
likely to introduce a bias in the measurements for high red- 
shift QSOs, in the sense of a steeper estimated continuum 
and a consequent underestimate of the Lyman absorption. Our 
low-redshift measurements, based on lUE+optical data, could 
therefore be more accurate but not homogeneous with respect 
to the high-redshift data, since the estimate of the "intrinsic 
continuum" is carried out on different rest-frame intervals. On 
the other hand, such a concern is probably less important than 
it may appear because a) in the present case, the only effect of 
the use of the optical+IUE data with respect to the lUE data 
only is to restrict the range of the possible continuum slopes 
to the steeper ones; b) the quasars used in the Da,b measure- 
ments have redshifts typically larger than those of the quasars 
used to construct the composite spectra, and for higher redshift 
QSOs the contribution of the Fe 2000 A complex is probably 
less important (e.g. Giallongo, Cristiani, Trevese 1992). 

To get an idea of the effects due to differences in the 
spectral ranges and/or resolutions it is interesting to note 
that the continuum slope estimated by us for PKS 0637-75 
(a = 0.30±0.05) is significantly flatter than the value obtained 
by O'Brien et al. (1988), 0.53±0.02, on the basis of lUE spectra 
and optical photometry, and that for MC1104+16 the spectral 
index obtained by Oke et al. (1984), 0.24 ±0.1, on the basis of 
very low resolution optical spectrophotometry, is significantly 
flatter than ours, 0.48 ±0.05. The above mentioned indices can 
also be compared with the "color-determined" values of Yu 
(1987) of 0.84 and 0.60 respectively. 

Other non-negligible sources of uncertainty are the match- 
ing between the optical and lUE flux calibrations and, espe- 
cially in the case of PKS 0637-75, the amount of galactic ex- 
tinction adopted. 

Since systematic errors are extremely difficult to evaluate, 
in Table 3 we have indicated as uncertainties in the q. Da, 
Db values only the "internal" uncertainties deriving from the 
different possible choices of the continuum slope and normal- 
ization. 



5. Comparison with models and conclusions 

We have compared the measured average depressions Da with 
the predictions derived on the basis of the model statistics of 
the Lye absorptions superimposed to a synthetic QSO spec- 
trum as in Giallongo, Gratton & Trevese (1990) and Giallongo 
& Cristiani (1990). 

As in the previous papers we have considered the following 
classes of absorbers. 

(i) The Lye lines not associated with metal line systems. 
These systems evolve in redshift with an average distribution 
of restframe equivalent widths according to the following equa- 
tion: 



d^N{z,W) 
dzdW 



#^<'-" 



exp(-W/W*) 



(2) 



A recent estimate of the Lye line distribution in the range 
1.8 < 2 < 3.8 and 0.2 <W < 0.5 A has been obtained from a 
small high resolution sample by Giallongo (1991). On average 
W* = 0.2 A and 7 = 2.1 was found. For W^ > 0.5 A a flatter 
distribution with W* = 0.3 A has been assumed following the 
equivalent width distributions obtained from large intermedi- 
ate resolution samples. A set of random absorption redshifts 
is extracted according to equation (2) with 7 = 2.1. We then 
associate to each redshift a random equivalent width extracted 
from a probability distribution P{W) oc exp(— H^/W^*) using 
the above W* values in the appropriate W ranges. 

(ii) hja lines associated with metal line systems. We adopt 
the same z, W distributions with 7 = 0.68 (Sargent, Steidel & 
Boksenberg 1989) and W = 1.05 A (Sargent et al. 1980) as 
in the previous papers. 

(iii) Lyman /3, 7, 5 are computed for each Lyo, using a 
standard curve of growth derived from a Voigt profile. We have 
adopted a Doppler parameter value of fo = 20 Km s^^ for weak 
lines {W < 0.5 A) and fe = 40 Km s"^ for strong fines {W > 0.5 

A). 

(iv) We have also included in the line statistics a popula- 
tion of damped systems with W* ~ 10 A to account for the 
average absorption observed at very high redshift (z > 4), (see 
Giallongo & Cristiani 1990 for details). For each emission red- 
shift we compute 100 synthetic QSO absorption spectra using 
Voigt profiles with natural damping only. 



Fig. 4. Comparison between the observations and the simulations for the average Da depressions. Crosses represent the average values 
from the Monte Carlo simulations. The contribution of the interstellar absorption lines has been subtracted to the Da values observed for 
PKS 0637-75 and MC 1104+16. 



The available data and simulations at various Zem for the 
average depressions Da are shown in Fig. 4. For Zem > 2.5 indi- 
vidual Da values are computed from optical spectra of QSOs 
taken by various authors at different resolutions (6-100 A). 
For Zsm < 2 average depressions are derived from lUE spectra. 

We have represented both the data of O'Brien, Gondhalekar 
& Wilson (1988) and Bechtold et al. (1984). We note that the 
values derived by Bechtold et al. are systematically higher, 
probably because of the low resolution adopted in the opti- 
cal spectra (80 A) which could cause an overestimate of the 
continuum level at wavelengths longer than the Lyo, emission. 

When comparing observed and predicted Da,b values on 
the basis of the Lja statistics, we have to remove the contribu- 
tion of the interstellar absorption lines, that is non-negligible 
when the average depressions are very low (e.g. < 0.1). For 
example, in the high resolution spectrum of 3C 273 (Bahcall 
et al. 1991), we can evaluate a value of Da — Wtot/AA = 0.02 
due to the interstellar absorption lines. We can estimate the 
analogous contribution for PKS 0637-75 and MC 1104-1-16 by 
considering the interstellar absorptions observed in the 3C 273 
spectrum in the range corresponding to the A region of PKS 
0637-75 and MC 1104+16, and scaling their equivalent widths 



according to the ratio between the EW of the 2798 A Mg II 
absorption feature observed in our quasars and the one in 3C 
273. In this way an extra Da value of 0.009 is derived for PKS 
0637-75 and of 0.006 for MC 1104+16. 

After this correction we have, at 2 ~ 0.65, an average ob- 
served value for the two QSOs, Da = 0.047 ± 0.013 which is 
larger than but still consistent, within the relatively large un- 
certainties, with the one predicted by the redshift evolution of 
the Lyo, lines Da = 0.026 ± 0.015. In this way, marginal ev- 
idence is found for a departure from the standard Lyman a 
evolution down to 2; ~ 0.65. The excess of absorption could be 
due to the presence of the diffuse component of neutral hydro- 
gen (the Gunn-Peterson effect) or to a slower evolution towards 
low redshifts of the Lya lines. 

Indeed recent HST high resolution spectra obtained at 
lower redshifts 2: < 0.3 seem to imply a slower evolution with 
7 ~ 0.7 (Bahcall et al. 1992). Adopting this evolutionary rate 
from 2 ~ 2.5 down to 2 = 0.65 we obtain Da = 0.047± 0.015 in 
agreement with our observed absorptions. Thus no significant 
Gunn-Peterson effect is detectable in the interval 0.65 < 2 < 
4.2. 



The implications of the change in the Une statistics would 
be important for the physical models of the Lye clouds and for 
the cosmological evolution of the intergalactic medium and the 
UV ionizing flux. In the standard scenario of pressure confined 
cloud model, the change in the cosmological evolution of the 
clouds is related to a quick drop of the UV ionizing flux at very 
low redshifts (Ikeuchi & Turner 1991). 

The uncertainties involved in the process of estimating the 
continuum slope and normalization, matching the spectra ob- 
tained in different spectral ranges with different instruments 
and detectors, evaluating correctly the galactic extinction are 
such to indicate that an improvement of the present estimates 
at low redshifts with low-resolution spectrophotometry will be 
extremely difficult. Higher resolution line statistics at low z is 
obviously the way to follow. 
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